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ABSTRACT 


Presented herein is an analog computer program capa- 
ble of predicting an orbiting spacecrafts external surface 
temperatures. By dividing a spacecrafts external surface 
into sections and applying transient heat transfer and ther- 
modynamic equations to the boundaries of each section, an 
analog circuit describing the total incident energy as a func- 
tion of spacecraft orbital position is written. The source of 
the external incident energy simulated can be from any com- 
bination of solar radiation, earth albedo, earth thermal 
emission, solar paddle emission and reflection. The sources 
of internal energy incident to the inner surfaces can be from 
instrument power dissipation, and from radiation inter- 
change between adjacent and opposite surfaces. Using the 
AOSO proposed spacecraft as the thermal model, the analog 
computer solution is shown for the parameters involved. 
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USE OF ANALOG COMPUTATION IN PREDICTING 
DYNAMIC TEMPERATURE EXCURSIONS 
OF ORBITING SPACECRAFT 


by 

Frank J. Cepollina 
Goddard Space Flight Center 


INTRODUCTION 


An analog computer program suitable for the analysis and design of spacecraft thermal 
control systems is described. The program was written essentially for the Advanced Orbiting 
Solar Observatory (AOSO) configuration; however, it can easily be adjusted to handle a multitude 
of configurations and orbital conditions. Given a set of thermal properties and orbital ther- 
mal flux, the program will yield skin temperatures as a function of the spacecraft orbital posi- 
tion. If optimum spacecraft surface thermal properties are desired so as to minimize temperature 
excursions, the program will yield the optimum thermal properties to achieve minimum tempera- 
ture excursions. In addition, the program is capable of solving internal temperature gradients 
across instrument packages, structural members, and thermal insulation materials as a function 
of spacecraft orbital position. Thus, detrimental temperature gradients resulting from various 
modes of operation (occulted, standby, failure) and various patterns of energy dissipation (on-off 
operation of instruments, batteries, and controls) can be predicted, and optimization of the thermal 
design achieved through the use of this program. 


THE PROGRAM MODEL 


The basic capabilities of the differential analog computer are the solution of ordinary 
differential equations. A typical spacecraft thermal analysis problem falls within the scope of 
these capabilities. Taking the AOSO spacecraft (4-solar paddle configuration) as an example 
(Figure 1), the cylindrical outer surface is sectioned by system boundaries, as shown in Figure 2. 
Each section of surface receives external environment thermal flux and internal dissipated energy, 
while simultaneously exchanging energy through reflection, emission, and conduction with space, 
spacecraft appendages, and the remaining sections of surface. The general differential energy 
balance equation, which can be written for each section, follows. 
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Figure 1— Three-axis stabilized cylindrical 
spacecraft (AOSO). 




^ solar 



Figure 2— Thermal model of a cylindrical spacecraft. 
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Earth emission flux 
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function of orbit position 
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where 


W = Mass of the cylindrical section (thermal shield section), 

C p = Specific heat of section material, 

Q ap = Solar energy entering through lens apertures on solar normal face, 
e s 3 External surface emissivity of thermal shield, 

F e s Earth to spacecraft thermal emission form factor— a function of spacecraft orbit 
and location in orbit ( f p ) , 

A s = Surface area of each section of the thermal shield, 

E t = Earth thermal emission (average) = 66.38 BTU/hr-ft 2 , 
a s h External surface solar absorptivity of thermal shield, 

F b Earth to spacecraft form factor for earth albedo— a function of spacecraft orbit and 
position in orbit with respect to the earth (f ), 

A p s Projected area of each section of thermal shield, 

S = Solar constant = 443 BTU/hr-ft 2 , 

a = Stefan Boltzmann's natural constant = 0.178 x 10' 8 BTU/hr-ft 2 (°R) 4 , 


P sp = Solar reflectivity of solar paddles, 

F v = Geometric form factor between solar paddles and each thermal shield section, 

(F pRs ) gr - Gray body configuration factor for thermal emission between paddle and thermal 
shield section, 


T 


P s 


T 

s 


T 

space 


^cond 


q 


rad 


in t 


Temperature of paddle surface, 

Temperature of thermal shield section: for each section T , , T T T * t, • 
t 2 , T 3 , T 4 respectively, 

Temperature of space, assumed equal to zero °R, 

Thermal conduction to or from each thermal shield section, 

Net emission exchange from the internal surface of each thermal shield section. 
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ASSUMPTIONS 


The analysis presented in this report for purposes of illustration is based on the AOSO con- 
figuration and its parameters, including: 

1 . The AOSO cylindrical spacecraft, 4-paddle configuration. 

2. The AOSO orbit (300 nm near-polar orbit): 

a. Continuous solar pointing spacecraft 

b. Space-oriented three-axis stabilized spacecraft. 

3. Anal og approximation of the earth albedo and IR flux curves as computed by the digital 
orbital mechanics program. 

4. Other assumptions as noted throughout the report. 


MECHANIZATION OF THERMAL MODEL 


For each thermal system boundary section of Figure 2, the general differential equation is 
rewritten in terms of coefficients and variables (References 1 and 2) and with the sign convention 
that energy entering a system boundary is negative and that leaving is positive.* 


For Section 1, 


-«.C P ^= -A 1 -B 1 -C 1 (f p )-D 1 (f p )-E 1 -F l (^ psl -'lt 1 )+G 1 (T^ 1 ) + I qcond(i3)(i4) + 2 q rad(i , 3 4 } • (1) 
For Sections 2, 3, and 4, the following equation applies, 

V ’ dW,.j.4 


^2 ,3,4 ®2,3,4 ^2,3,4 ®2,3,4^p) ^2,3,4 


^2,3,4 02. 2,3,4 ^s2,3,4^ + ^2,3,4 v ‘s2,3,4 


(^s2.3,4^ + ^ t lcond 2 3 4 + 3 4 


( 2 ) 


(f p ) 

('wJDO 

rs t 

+ T S 1(F, +G,) 


-(A, +B, + E, ) 


T P$I ( F 1 ) 


{ - D ><V 



3 1 ) 

( 



SUMMATION OF ALL 
EQUATION CON- 
STANTS TO =WCp dl 
dr 


INTEGRATING dT S i/dr 
TO GENERATE T S] AS 
A FUNCTION OF TIME 


DOUBLE SQUARING 
T SI TO GENERATE T^ 


Figure 3 — Basic circuit for each section. 


The constant coefficients of these equations 

are A x 2 3 4 ; ^ 234 * ^ 1 , 2 , 3 , 4 ’ ^ 1 , 2 , 3 , 4 ? 

G 1 2 3 4 . The variable coefficients are 
c i! 2 ,3,4 ( f P ) and D i, 2,3,4 ( f P )- Both of the 
variable coefficients are dependent on the 
specific spacecraft orbit and the position of the 
spacecraft in that orbit (relative to the earth) 
at any instant of time. Each of the above 
equations can therefore be mechanized by the 
following differential analog circuit (Figure 3). 

By using four analog circuits identical to 
that in Figure 3, (one for each of the four sys- 
tem boundaries equations), it is possible to 


*This sign inversion facilitates efficient analog mechanization of the general equations and insures solution stability as per Rouths 
criterion (Reference 3). 
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generate the temperatures of each of the four sections as a function of time. However, the result- 
ing temperatures would be generated independently of any conduction or internal radiation exchange 
between sections and internal components of the spacecraft. In reality, spacecraft do interchange 
energy across these theoretical system boundaries via conduction and internal radiation. There- 
fore, it is essential that a means of conduction and radiation coupling between each of the four cir- 
cuits be included in the program. 


Conductive Coupling 

The net thermal energy conducted to any section is a function of the temperature difference 
between any given section and its adjacent sections. 

For example, for Section 1, if the thermal conduction resistance, R, is assumed equal for each 
of the four system boundaries, then 


In 

M co 


- — (T. 
R V 4 


V +^(T t -T a ) 


2 

R 



(3-1) 


Similarly, for Sections 2, 3, and 4, 




— T, - 1- T, 
R 3 R 4 


(3-2) 


2 


q 


con d 


3 




(3-3) 


2 



(3-4) 


By incorporating the above equations into 
four circuits (each identical to the circuit 
described in Figure 3), conductive coupling of 
each section is accomplished. The mechani- 
zation of the equation is achieved through the 
use of feedback loops (Reference 3), and the 
resulting interconnected circuits are described 
in Figure 4. 


Internal Radiation Coupling 

The technique used to incorporate internal 
radiation coupling into the program is very 


{ A - 1 ) 


(fw) — Q 

— Q ( _ d 


(A-2) 


(A-3) 


(A-4) 



— O 


i-i^ 


(A-2) 


Figure A — Basic circuit diagram 
with conduction coupling. 
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much the same as that used to achieve con- 
ductive coupling. A simplified internal config- 
uration of the spacecraft experiment compart- 
ment is used to demonstrate the technique. 
Figure 5 illustrates the thermal shield (with 
the four arbitrary system boundaries and a 
hexagonal central support column or tube). In 
addition, a solar normal cover and an aft cover 
are added to the thermal model. Both forward 
and aft covers are assumed to be at constant 
temperatures, and conductively insulated from 
the thermal shield. These assumptions reduce 
the complexity of the program, but are not 
requirements to the operation of the program. 

The net radiation exchange to Section 1 from 
all other sections, central tube, and forward 
and aft covers, can be expressed by the equation 


Figure 5— Internal thermal model. 


^radj = Fi4 A ^(r> - TJ) 4 F,, Aa(T? - T«) 
+ F 12 Aa(T? -TJ) + F 15 Aa(Tj -T«) 
+ F 17 Aa(T? -T}) +F 16 A CT (T?-TS) 
+ F 16 Aa(T? -T? a ). 


(4) 


For the particular configuration, F 12 is approximately equal to zero and by symmetry, f 13 = f 14 
and F 1S = F 16 . Combining and expanding terms yields the following equation: 


2q radl =F 14 A 1 a(2T< - TJ - Tf) 

+ F 15 A 1 <r(2Tt - T^) 

+ F l7 A i-(Tt -lt A ) . 

By analogy, each of the remaining section equations are written 

ad 2 =F 2 4 A 2 ^( 2 ^-' l 4 4 -^) 

+ F 25 A 2 a(2T^-lt-lt) 

+ F 27 A 2 cr (^2 “ ’ 


(4-1) 


(4-2) 
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S ^ad 3 = F 31 A 3-( 21 1-' l4 2- T t> 

\ 

+ F 35 A 3°-< 2T ?- 1 1- T 1) 

+ F 37 A 3 CT ^- T 7C)- (4-3) 

and 

4 

+ F 45 A 4^C2Tt -T< -T«) 

+ f 47 a 4 o- ( tJ - t^ d ) . (4-4) 


To implement placing these equations into 
the circuits of Figure 4, the technique of feed- 
back loops is again used. However, at this point 
no circuit exists to generate T 7A , T ?b , t ? c , t 7D , 
T. andT*. If it is assumed that these values 
are known constants, then mechanization can 
proceed. If they are known independent vari- 
ables, then a function generation of the inde- 
pendent variables is required as inputs to the 
circuits of Figure 4. On the other hand, if they 
are unknown dependent variables (functions of 
the temperature and energy received from the 
surroundings), then it is necessary to write 
additional differential equations and mechanize 
circuits to generate these variables. The tech- 
nique used to write and mechanize the addi- 
tional equations is identical to that described 
in the preceding pages. 


(A- 1 ) 

( A- 4 ) - 
( A- 3 ) - 
(A-l) ■ 
CONST. 
CONST. 

(A-l) - 
CONST. 

( A- 2 ) 
(A-4) ■ 
(A-3) ■ 
(A-2) ■ 
CONST. 
CONST. 

(A-2) - 
CONST. 


(A-3) 
(A-2) 
(A-l ) 
(A-3) 
CONST. 
CONST. 

(A-3) 

CONST. 



(B-l) 


(B-2) 


(B-3) 


So as not to confuse the illustration of 
the analog computation technique, the terms 
T 7 a b c d y T s T e be considered con- 

stant, even though it is not a requirement to 
the analog solution. 

Figure 6 illustrates the mechanization of 
these radiation coupling equations, 4-1, 4-2, 



Figure 6 — Mechanization of 
radiation coupling equations. 
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4-3, and 4-4. Note that for circuit B-l, T*, T* , 
and T* , functions are generated by circuits 
A-l, A-4, and A-3, respectively. Multiplica- 
tion by factors of two of the T^| function is 
achieved by using gains of two at each respec- 
tive summing amplifier input. By assumption, 
T^, T*, and T^ a are constants. Now, by taking 
the output of the final summing amplifier and 
inputting this into the summing amplifier of 
circuit A-l, the radiation exchange (coupling) 
loop is completed for Section 1. By the identi- 
cal procedure, radiation coupling is accom- 
plished for Sections 2, 3, and 4. Figure 7 
illustrates the augmentation of radiation coup- 
ling into the main circuits. 


SIMULATION OF TIME 
VARYING HEAT FLUXES 


The generation of earth thermal emission 
and earth albedo heat flux is accomplished by 
simulation rather than by solution of the 
orbital mechanics equations. Although it is 
possible to solve these equations by analog 
techniques, the higher degree of accuracy of 
digital computation techniques cannot be disputed. Therefore, this particular technique simulates 
the time varying heat flux as predicted from the digital orbital mechanics program. 

The approach taken to develop the simulation is as follows (Reference 4): 

1. Write the Fourier Series expression which approximates the shape of the albedo and IR 
flux curves for each of the four boundary sections. 

2. Compute the amplitude and phase shift required for each flux curve for each of the four 
sections. 

3. Mechanize a single cosine/sine generation circuit. 

4. Mechanize each of the Fourier expressions through the use of diode chopped feedback 
amplifiers and the single sine/cosine circuit. 

5. The potentiometer settings are then adjusted so that the output of each amplifier matches 
its respective digital flux curve. For the ease of making rapid changes in the external emissivity, 
e, and solar absorbtivity, a, of each spacecraft section, a separate row of potentiometers is placed 
in the circuit. 



Figure 7— Circuit diagram incorporating both 
radiation and conduction coupling. 
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Figure 8 illustrates two flux generation 
circuits coupled to a single sine/cosine gen- 
eration circuit. Figures Al, A2, A3, and A4 of 
Appendix A present the simulated earth emis- 
sion and albedo fluxes for each of the four 
sections of spacecraft surface. To demon- 
strate the accuracy of the simulation, a 
comparison is made between computed digital 
values of flux (Reference 5) and the values 
simulated by the flux generation circuits (for 
each of the four surface sections). As can be 
seen in the referenced figures, the simulation 
can be made to match the digital solution 
closely. The circuits which generate these 
time varying fluxes are shown in Figure A5 of 
Appendix A. 

TEMPERATURE ANALYSIS (COMPUTER RUNS) 


The circuit diagram of Figure A5, Ap- 
pendix A, represents the complete mechani- 
zation of the preceding equations into a dif- 
ferential analog computer program. Through 
the use of this circuit, the program is demon- 
strated. The particular temperature analysis 
(computer runs) is made for a cylindrical, 
48-inch diameter, solar pointing, three-axis 
stabilized spacecraft in a 300 nm near-polar 
orbit. The sample evaluation of equation 
constants, coefficients, and scaling is included 
in Appendix A. The earth flux condition is 
taken as the summer solstice (June 21) orbital 
period. The simulated flux for each node is 
described in Appendix A, Figures Al, A2, A3, 
and A4. 

The following figures illustrate computer 
solutions for various enumerated thermal 
properties of the spacecraft’s surface (Fig- 
ures 9, 10, 11, and 12). 



Figure 8— Sample earth thermal 
flux simulation circuits. 
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90 f— 


3 

H- 

2 
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48-INCH DIAMETER CYLINDRICAL VEHICLE; 10 FT IN LENGTH 
EXTERNAL EMISSIVITY, e =0.08 
SOLAR ABSORPTIVITY, a =0.268 

SKIN CONDUCTION, 1/R = 0.5117 BTU / hr - °F/ section 
INTERNAL RADIATION BLOCKAGE = 100% 

INTERNAL POWER DISSIPATION, Q = 79.5 BTU /hr/ section 
ONE ORBIT = 96 MINUTES 



TIME (minutes) 


Figure 9— Analog computer run for 
thermal design optimization. 


9 





100r 


CONDITIONS: 


l 


90 1 — 48- INCH DIAMETER CYLINDRICAL VEHICLE 

10 FT IN LENGTH 
EXTERNAL EM1SS1VITY, e =0.08 
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80 1— SKIN CONDUCTION, l/R * 2.0468 BTU/hr - °F/secHon 
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INTERNAL POWER DISSIPATION, Q = 79.5 BTU/hr /section 
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Figure 10— Analog computer run for thermal 
design optimization. 


Figure 11— Analog computer run for thermal 
design optimization. 


CONCLUSION 
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The basic capabilities of the analog com- 
puter comprise the solution of ordinary dif- 
ferential equations. A typical spacecraft 
thermal analysis problem falls within the 
scope of these capabilities. Using the proposed 
Advanced Solar Observatory as an example, 
thermal analysis solutions have been demon- 
strated. The thermal model consisted of the 
spacecraft’s cylindrical surface sectioned by 
systems boundaries, and the general differen- 
tial energy balance equations were written for 
each boundary section. It has been further 
demonstrated that a differential analog circuit 
can be designed to generate the solution of 

each differential energy balance equation. Time varying orbital flux can also be accurately simu- 
lated, and the flux simulating circuits were designed so that thermal properties (surface emissivity, 
absorptivity, form factors) can be readily varied to accommodate rapid optimization of thermal 
properties. Finally, several analog runs were made to demonstrate the capabilities of the 
program. 


60 80 
TIME (minutes) 

Figure 12— Analog computer run for thermal 
design optimization. 


Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland, May 2, 1966 
829-11-75-01-51 


10 




































REFERENCES 


1. Kreith, F., "Radiation Heat Transfer for Spacecraft and Solar Power Plant Design," New York: 
International Textbook Company, 1962. 

2. Powers, E. I., "Thermal Radiation to a Flat Surface Rotating About an Arbitrary Axis in an 
Elliptical Earth Orbit: Application to Spin Stabilized Satellites," NASA Technical Note 
D-2147, April 1964. 

3. Rodgers, A. E., and Connoly, T. W., "Analog Computation in Engineering Design," New York: 
McGraw Hill Book Company, 1960. 

4. Jackson, A. S., "Analog Computation," New York: McGraw Hill Book Company, 1960. 

5. Stevenson, J. A., and Grafton, J. C., "Radiation Heat Transfer Analysis for Space Vehicles," 
Aeronautical Systems Division Technical Report 61-119, U. S. Air Force, December 1961. 


11 




Appendix A 


Sample Preparation of Input Data for 
Spacecraft Temperature Analysis by Differential 
Analog Computer 


General Equation 

For a three-axis stabilized spacecraft, cylindrical in configuration, solar pointing in orienta- 
tion, and in a 300 nm near-polar orbit, the following general equation applies: 

( WC P dr )i, 2> 3 I N _A 1.2,3,N +B 1,2,3,N + C 1 , 2 , 3 , N ( f p ) + , 2 , 3 , N ^ +E 1,2,3,N +F l,2,3,N^ 1 ps 1>2>3 ,N “ I. 2, 3,N^ 


G 1 , 2 , 3 , N ( T s 1>2t3 ,4 ) f ^ Q<:ond 1 ,2.3. 4 ± ^ Qrad int l,2,3,N * (Al) 


The subscripts 1, 2, 3, and N represent the number of respective equations for N thermal boundary 
sections (nodes) comprising a spacecraft thermal model. If the spacecrafts shell is cylindrical* 
and is divided into four thermal boundary sections, then four equations identical to Equation Al 
will describe the temperature history of the spacecraft’s shell. Since each equation is identical 
for the cylindrical thermal model being studied, only the evaluation of one equation describing one 
section (node) will be discussed herein (the technique being the same for each remaining equation). 


Evaluation of Terms 

For the aforementioned three-axis stabilized spacecraft, the constant terms of Equation Al are 

A = Q^, solar energy absorbed through apertures, 

B = q int , internal electrical power dissipation, 


and 


E = Solar paddle reflections absorbed by skin section. 

The equation independent variable terms are 

C(f ) = Incident earth IR flux (simulated functions of sin and cos of e ), 

and 

D(f p ) = Incident earth albedo flux (simulated functions of sin and cos of 0 ). 


13 



The equation dependent variable terms are 
F(T* s -T*) = Paddle-to-skin IR exchange, 

G(T 4 ) s Skin surface emission to space, 

I q cond s Skin section conduction, 
la = Skin section internal radiation exchange, 

^rad int 

and 

dT/dr = Transient change in temperature of a section with respect to time. 

The computet variables of Equation Al, Section 1, are 

Equation Variables Computer Variables 


QA S1 

dr 


T.i 


T T , 

A sl’ A s3 

T S 4 

T .i- 

T s3’ T s4 

j4 

1 sl* A s 3 ’ 

^4 

T 4 

A sl’ 

T 4 , T 4 

A s3 A s4 

6 


cot 



Computer Equation 

Rewriting Equation Al in terms of the computer variables and using the sign convention that 
energy entering a system boundary is negative and that leaving is positive, yields 

- WC p [T sl ] = - Aj - Bj - Ej - Cj (f p [cos cat] ) - Dj (f p [cos wt] ) + (F+G) [T^] 


-F05.,) + -^—(2 [T sl ] - [T s4 ] - [TJ) 

■**cond 

+ (2 W ~ tTl 4 ] - [T^ ) 


+ (2 [tjj > - -i- cn - 

R rad 2 K rad 2 


+ 


1 





(A2) 
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Inputting values from Table 1 into Equation A2 yields the following unsealed computer equation: 


4. 8 [T s ] - - 267.33 -63-44 { f [cos o;t ] } - 48.33 { f [cos cot] } 


+ 0.3741 x 10' 8 { [T? ] } + 0.5117 (2[T ] - [T, ] - [T a ] } 

S 1 S 1 b 4 s 3 


+ 0.3988 x 10‘ 8 {2 [T 4 ] - [T 4 ] - [T 4 ] } 

S 1 s 4 s 3 


+ 0.3551 x 10‘ 8 {2 [T 4 ^ - 2(5 2 0) 4 ) + 0.8915 x 10‘ 8 { [T 4 ^ -(520) 4 } 

The amplitude scaled equation is 


(A3) 


0.48[T sl ] 


U-i 


i 

l 

LV io/ j 

J 


+ 0.05117 


s i 

10 


10 


s 3 

10 


0. 03998 < (2 X 10“ 4 ) 


+ 0.03551 <(2 x 10' 4 ) 


fefl - 

/T \ 4 

I 111] 

\ 10 / 


(10- 4 ) 


s 4 

10 / 


) -( 10 - 4 ) 


3?) 


51-8 




+ 0.08915|(10- 4 ) 


-65.0 


(A4) 


Time scaling the computer integration and earth flux generation rates, 


= 0.48 

10 


dr 


lo.«il , 

10 0.48 


_L = 0.20833 [0.48T ] °F/sec , 
10 
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Table 1 


Computer Program Input Information. 


Problem Variable 

Computer Variable 

Scaled Variable 

dT Sl 

dr 

T 

[0.48T] 

V V X *4 

T s * T s T s * 

s 1 s 3 s 4 

ft] 

T s 3 ’ 

X 4 t 4 * T 4 * 

s l’ s 3 * s 4 

[■”-'(¥] 

e 

cut 

[cat] 

Equation Constants 
and Coefficients 

Problem Constants 

Value 

A, 

Qap 

15.38 BTU/hr 

B i 

^int 

64.12 BTU/hr 

E ! 

p a F A S 

' sp s v p 

48.33 BTU/hr 

F i( T p s ) 

(F ps ) gr A s O- (T ps ) 

139.50 BTU/hr 

C l 

£ s F A E * 

63.44 BTU/hr 

D ! 

a s F pA p S 

48.33 BTU/hr 

F l 

(FpAr A s * 

0.1557 x 10' 8 BTU/hr(°R) 4 

G i 

€ s A s CT 

0.2184 x 10' 8 BTU/hr(°R) 4 

(F + G)j 

CT t(F ps ) gr A s + £ s AJ 

0.3741 x 10' 8 BTU /hr(°R) 4 

1 

AX 

0 SI 17 BTU /hr °F 

^cond 

KA 


1 

ft rad int 

F 14 A i a ’ F lS A l CT - F 17 A l°' 

0.3988 x 10" 8 , 0.3551 x 10‘ 8 , 0.8915 

wc p 

WC p 

4-8 BTU/°F 

T 5 

T 

1 s 5 

5.20 °R 

T 6 

T 

s6 

5.20 °R 

T 7A 

T 

1 7 A 

5.20 °R 


'Generated by Equation Al, written for thermal boundary sections 3 and 4. 



and 


T 

— = 12.500 [0.48T] °F/min. = computer integration rate. 

10 

To be compatible with the computer integration rate cot ~ 0.48 radians/orbit revolution, 

t = 1.6 hrs. per orbit = 96 minutes/orbit. 

Therefore, [oil = 0.0655 radians/ min, flux generation rate. 

Scaled Analog Diagram 

In a technique identical to that above, three additional equations are evaluated, scaled, and 
mechanized. Table 2 lists the scaled coefficients for all four equations in terms of potentiometer 
and gain settings for the completely mechanized analog circuit (Figure A5). Figures Al, A2, A3, 
and A4 represent the earth flux generated for each node by the flux simulation circuits. 


17 



Table 2 


Potentiometer Assignment Sheet. 


DATE: SAMPLE PROBLEM 


P00-P59 

PROBLEM: THERMAL ANALYSIS CIRCUITSSIS-1096 


POT. 

NO. 

P 

PARAMETER 

DESCRIPTION 

SETTING 

STATIC 

CHECK 

SETTING 

RUN 

NUMBER 1 

POT. 

NO. 

P 

PARAMETER 

DESCRIPTION 

SETTING 

STATIC 

CHECK 

SETTING 

RUN 

NUMBER 1 

00 

Solar Face Temp. 

.5200 

.5200 

30 

Scaling (gain = 12.5) 

.5000 

.5000 

01 

Aft Cover Temp. 

.5200 

.5200 

31 

Initial Temp. -Sur. 4 

.5000 

.5000 

02 

Tube Temp. 

.5200 

.5200 

32 




03 

Plotter 


.5000 

33 

Conduction-Skin 

.1022 

.1022 

04 

IR Flux f(F e ) Control 3 

.5249 

.5249 

34 




05 

Plotter-X Drive 


.0096 

35 

Initial Temp. -Sur. 3 

.5000 

.5000 

06 

Sur - 1 9 rod in, 

.03998 

.03998 

36 

Q Internal Sur. 3 

.2673 

.2673 

07 

Sur ' 1 9r=d int 

.03551 

.03551 

37 

A s Flux-Control-Sur. 3 

.0595 

.0595 

08 

Sur ' 1 9 rad int 

.08915 

.08915 

38 

A s Flux-Sur. 3 (a s ) 

.2680 

.2680 

09 




39 

Conduction-Skin 3 

.2044 

.2044 

10 

Scaling (gain = 12.5) 

.5000 

.5000 

40 




11 

Initial Temp. -Sur. 1 

.5000 

.5000 

41 




12 

A s Flux Control-Sur. 1 

.0083 

.0083 

42 




13 




43 




14 

IR Flux Control-Sur. 1 

.0900 

.0900 

44 




15 

Sur - 2 9 rod int 

.03998 

.03998 

45 




16 

Initial Condition 

.1000 

.1000 

1 46 

Conduction-Skin 

.1022 

.1022 

17 

Sur - 2 9 rod int 

.03551 

.03551 

47 




18 

IR Flux Control-Sur. 3 

.0456 

.0456 

48 

Conduction-Skin 

.1022 

.1022 

19 

IR Flux f(F e ) Control 2 

.8830 

.8830 

49 




20 

Scaling (gain = 12.5) 

.5000 

.5000 

50 




21 

Initial Temp. -Sur. 2 

.5000 

.5000 

51 




22 

A s Flux-Sur. 2 (ci s ) 

.2680 

.2680 

~52 




23 

IR Flux Control-Sur. 2 

4)864 

.0864 

_ 53 




24 

A s Flux Control-Sur. 2 

.0126 

.0126 

54 




25 

IR Flux Control-Sur. 4 

.0456 

.0456 

55 




26 

IR Flux f(F e ) Control 4 

.5249 

.5249 

56 




27 

A s Flux Control-Sur. 4 

.0456 

.0456 





28 

A s Flux Control-Sur. 4 

.7225 

.7225 

58 




29 




59 
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Table 2 (Continued) 


Potentiometer Assignment Sheet. 


Q00-Q59 

DATE: SAMPLE PROBLEM PROBLEM: THERMAL ANALYSIS CIRCUIT SSIS-1096 


POT. 

NO. 

Q 

PARAMETER 

DESCRIPTION 


SETTING 

STATIC 

CHECK 

SETTING 

RUN 

NUMBER 1 

POT. 

NO. 

Q 

PARAMETER 

DESCRIPTION 

SETTING 

STATIC 

CHECK 

SETTING 

RUN 

NUMBER 1 

00 

Plotter 



.5000 

30 




01 

Plotter 



.5000 

31 

Conduction-Skin 4 

.2044 

.2044 

02 

Plotter 



.5000 

32 

Q Internal-Sur. 4 

.2673 

.2673 

03 

Plotter 



.5000 

33 

IR Emission-Sur. 4 

.3742 

.3742 

04 

IR Flux f(F e ) Control 

1 

.4225 

.4225 

34 

Plotter 


.5000 

05 

A s Flux f(F p ) Control 

1 

• 2720 

.2720 

35 

A s Flux f(F p ) Control 2 

.2567 

.2567 

06 

Sur ' 3 Irod in. 


.03998 

.03998 

36 

IR Emission-Surf. 3 

.3742 

.3742 

07 

Sur - 3 1 rod in. 


.03551 

.03551 

37 

Scaling (gain =12.5) 

.5000 

.5000 

08 

Sur - 3 <1 rod in. 


.08915 

.08915 

38 

A s Flux Control-Sur. 3 

.0101 

.0101 

09 

A s Flux-Sur. 1 (a s ) 


.2680 

.2680 

39 

A 5 Flux f(F p ) Control 3 

.2605 

.2605 

10 

IR Flux-Sur. 1(d) 


.0800 

.0800 

40 

Ave. Temp. Divider 


.3000 

11 

Conduct Ion -Skin 1 


.2044 

.2044 

41 




12 

Q Internal-Sur. 1 


.2673 

.2673 

42 




13 

IR Emission-Sur. 1 


.3741 

.3741 

43 




14 

A s Flux Control Sur. 

1 

.7700 

.7700 

44 




15 

Sin & Cos Generator 



.0655 

45 




16 

Sin & Cos Generator 



.0655 

46 

Conducti on 

.1022 

.1022 

17 





47 

Sur - 4( i rod int 

.03998 

.03998 

18 

Sur - 2 ^rad in. 


.08915 

.08915 

48 

Sur - 4 9 rod in. 

.03551 

.03551 

19 

IR Flux-Sur. 3(d) 


.0800 

.0800 

49 

Sur - 4 9 rod in. 

.08915 

.08915 

20 

IR Flux-Sur. 2(d) 


.0800 

.0800 

50 




21 

Conduction-Skin 2 


.2044 

.2044 

51 




22 

Q Internal-Sur. 2 


.2673 

.2673 

52 




23 

IR Emission-Sur. 2 


.3741 

.3741 

53 




24 

A s Flux Control Sur. 

2 

.7700 

.7700 

54 




25 

A s Flux-Sur. 4 (a s ) 

.2680 

.2680 

55 




26 

IR Flux-Sur. 4(d) 


.0800 

.0800 

56 




27 

A s Flux f(F p ) Control 

4 

.2630 

.2630 

57 




28 

A s Flux Control-Sur. 

4 

.9100 

.9100 

58 




29 





59 
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Figure Al — Incident thermal flux to 
spacecraft surface section 1. 
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Figure A2 — Incident thermal flux to 
spacecraft surface section 2. 
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Figure A3— Incident thermal flux to 
spacecraft surface section 3. 
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Figure A A — Incident thermal flux to 
spacecraft surface section 4. 
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